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Abstract 

New exact solutions of the steady and incompressible 2D MHD equations in polar 
coordinates are presented. The solutions describe the process of reconnective mag- 
netic annihilation in a curved current layer. They are particularly interesting for 
modeling magnetic reconnection in solar flares caused by the interaction of three 
photospheric magnetic sources. 
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1 Introduction 

Magnetic reconnection is a fundamental process in many areas of plasma 
physics. In particular it is proved to play a key role in active phenomena 
such as solar flares and geomagnetic substorms [1,2]. Since the pioneering 
works of Parker [3], Sweet [4] and Petschek [5], several exact solutions describ- 
ing reconnective annihilation in two dimensions in Cartesian coordinates were 
discovered e.g. by Sonnerup and Priest [6], Craig and Henton [7] and Priest 
et al. [8]. The aim of this paper is to present similar solutions for a curvilinear 
geometry in a form that makes it possible to model a class of solar flares. 
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2 Basic equations 



We consider stationary and incompressible plasma flows with uniform density 
and resistivity which are described by the equation of motion 

{v-V)v = -Vp+{VxB)xB (1) 

and by the Ohm's law 

E + vxB^r]WxB, (2) 
where the velocity v and magnetic field B satisfy 

(3) 

V-S = 0. 

The equations (1) and (2) are written in a dimensionless form. The distances 
have been rescaled with respect to a characteristic length Lg. The magnetic 
field and the plasma velocity have been normalized respectively to B^ and 
t'Ac which arc characteristic values of the magnetic field and of the Alfvcn 
velocity. The nondimensional resistivity rj corresponds to the inverse magnetic 
Reynolds number. The current density j is determined by the Ampere's law 

j^VxB. (4) 

Assume that the velocity and magnetic fields lie in one plane and do not 
depend on the z-coordinate perpendicular to this plane. Then the electric 

— * 

field E is uniform and parallel to the 2;-axis. In the plane where the magnetic 
and velocity field lie we use polar coordinates (r, 6) related to the Cartesian 
coordinates (x, y) as follows 

X — rsvaO, y — t cos 9 — d, (5) 

where d > so that the pole is below the line y — 0. The above mentioned 
Le is the distance from the origin of the coordinate system. Further we re- 
strict our consideration to the domain (y > 0,r < 1) in the corona, since the 
subphotospheric fiows are described by other equations which at least have to 
take into account the gravity. Due to (3) the velocity and magnetic fields can 
be expressed in terms of stream and flux functions ip and A, respectively. In 
polar coordinates this yields 

Eqs. (1) and (2) rewritten in terms ofip and A, have, respectively, the following 
form: 

[V',VV] = [AV2A], (7) 



2 



r 



Here the Poisson brackets are defined as 



^■''^^ drde drde' 



(8) 



(9) 



3 Form of the solutions 



We seek solutions in the form 

A{r,e) = A^{r)d + Ao{r), 



(10) 



^(r,^) =^i(r)^ + ^o(r), (11) 

where Ai, V'li ^0 and ipQ are unknown functions of r. Then the radial and 
azimuthal components of the magnetic and velocity fields are 



S,(r) = ^, Be{T,e)^-A^e-A,\ 

T 



(12) 



in which the symbol / indicates the derivative with respect to r. 
By substituting the ansatz (10) and (11) into (7) and (8) one obtains expres- 
sions depending linearly on 9. This means that in each of these expressions the 
coefficient of 9 and the sum of the terms not depending on 9 have to vanish. 
This yields a set of four ordinary differential equations which splits into two 
subsystems. The first of them is nonlinear and includes only Ai and -01, so 
that 



(rA/) - A 



(rA/) 



V'l>li-V'i>li+?7(A' + rAi") = 0. 
The second subsystem is linear in and ■00 and it has the form: 



— (r^)! ) - -01 
r 



-(rV'o')' 
r 



^{rA^)-A,\krA,') 
r ir 



(13) 
(14) 

(15) 



E + -[ij'^A^ - + ^(Ao' + rV)] = 0. 



(16) 



Since the obtained set consists of four equations for four unknowns, the as- 
sumed ansatz is compatible with equations (7) and (8). 
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4 Ideal solutions 



Let us consider first the case of ideal MHD, that is when 77 = 0. In this limit 
it is easy to see from eq. (14) that ipi must be proportional to Ai, that is 

A = (17) 

where a is an arbitrary constant. By using this result in (13) we obtain the 
equation 

Ai" + — ± A'Ai = 0, (18) 
r 

where A = — {Bj-'^dj /d9^ | Let us consider the case of a 7^ ±1. If A = 0, 
so that the current density does not depend on 9, the solutions of the system 
for vanishing resistivity are given by 

V^i = Cilnr + C2, = — (ci Inr + C2), (19) 

a 



, , a Er ah , , 1 Er h 

= r— ^ ^ + -r + — , ^0' = ^— ^—^ ^ar + -, 

CK^ — 1 Ci m r + C2 a ar o;^ — 1 Ci m r + C2 r 

(20) 

where Ci, C2, a and h are arbitrary constants. From (20), (12) and (4), it is 
easy to see that in the presence of a non-vanishing electric field E the 6- 
components of the magnetic and velocity fields as well as the current density 
have a singularity at the radius r = exp(— C2/C1) = rc, which is further on 
called the critical radius. In ideal MHD this is an indication of reconnection. 
The singularities disappear at E' = and in this case the arc r = turns into 
the separatrix line emanating from the null point (rc,0). Thus, the electric 
field is responsible for the appearance of the singularities at the separatrix: 
in our solution it drives the shearing flows across the other separatrix line, 
which is exactly the situation where the magnetic flux has to pile up at the 
separatrix aligned with the shear [9] . The flrst example of such a solution in 
a more simple geometry has been discovered by Craig and Henton [7], which 
showed also that the corresponding singularity is resolved by resistivity. It is 
shown below that this is valid for our solution as well. 



5 Resistive solutions 



To resolve the above singularity we adopt the method of matched asymptotic 
expansions by analogy with the work of Priest et al. [8], where it is used for 
solving a similar problem in Cartesian rather than cylindrical geometry. One 
can use such a method, because the dimensionless resistivity 77 is very small in 
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the solar corona as well as in many other astrophysical and laboratory plasmas. 
To solve the problem in this case the domain is separated in two different 
regions, a narrow layer enclosing the singularity and the rest of the plane. 
In each of the regions we first find the appropriate asymptotic expansions of 
the resistive solutions, then we match and combine them into a composite 
solution, which is approximately valid in both regions. 

The boundary conditions assumed for the unknowns in our problem are the 
following: 

Ai(l)=Sre, M^)=Vre, (21) 

^i(rc) = 0, Vi(rc) = 0, (22) 

Ao'in) = 0, (23) 
where B^^ and tVc are some values of radial components of magnetic and veloc- 
ity fields, respectively In the outer region (|r — rd > s/rj) the resistive terms 
in (14) and (16) are small and so the outer solution coincides in the lower or- 
der approximation with the ideal solution considered in the previous section. 
Using boundary conditions (21) and the above definition of Tc, eq. (19) can be 
written as 

= ^, = -^ln(^). (24) 

mTc Vrc/ luTc Vrc/ 

In this form, the outer solution (24) satisfies the inner boundary conditions 
(22) and so it can be used as an inner solution. A similar fact has been observed 
by Priest et al. [8] in the "Cartesian" analogue of our solution. The problem 
thus reduces to find the inner solutions for Aq and ipo . We assume that the 
following expansions 

Ao' = Al' + ^Al', V^o' = ^'o + Vvi^o, (25) 

are valid in the inner region: These expressions are inserted into eq. (15) 
and (16) where Ai and ipi are replaced by the first three terms of the series 
expansions of (24) about r^. The two resulting equations are then rewritten 
in terms of the inner variable 



r — Tc 




The terms with the same powers of r] are gathered in these equations to equate 
separately to zero the corresponding coefficients of the two lowest powers of 
77. This yields us four equations for the unknowns Aq', Aq', ifj^ and V'o'- Their 
solution, rewritten in terms of the variable r, determine approximate expres- 
sions for Aq and iJjq in the vicinity of Vf.. The latter is expanded then by small 
T] aX a, fixed r and matched with the series expansion of the outer solution 
about r = r^- Finally, the matched outer and inner solutions are combined 
into the following composite solution: 
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Fig. 1. Magnetic field lines (solid) and streamlines (dotted) for d = 0.05, E = 0.5, 
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CS, = 0.9, rc = 0.2, a = and 6 = 1. The plots of the field 



lines are superimposed to the distribution of the current density in gray half-tones. 
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where k — {vre^ — Bj-J^) / {vreTc^ ^{tc)) must be positive to avoid unphysical 
divergence of the solution at 77 — > 0. 

The corresponding magnetic and velocity field lines are plotted in Fig. 1. One 
can see from this plot that the separatrix line r = rc is the same for both 
sets of lines, while the other separatrices are different. The separatrices inter- 
sect at the magnetic null point (r = Tc, ^ = 0) and at the stagnation point 
(r = rc,0 = (arc + b/ r c) kr c^\n /v^-c'^). Such a structure implies the presence 
of a shearing component of the fiow parallel to the first separatrix r = Tc 
and transverse to the second magnetic separatrix. As mentioned above, this 
is a reason for the current layer formation along the first separatrix, which is 
confirmed by the corresponding distribution of current density in Fig. 1. 
In comparison with its "Cartesian" analogue [8] the obtained solution is much 
less symmetric. It describes the plasma fiow in curvilinear magnetic configura- 
tion with an arc-like current layer separating a dipole-like structure from the 
surrounding unipolar field. This is particularly interesting for modeling mag- 
netic reconnection in solar flares. In fact there arc observational evidences that 
many flares occur in conflgurations with three magnetic flux concentrations 
on the photosphere [10]. The proof of this conjecture as well as the detailed 
investigations of other solutions with A 7^ will be presented in a forthcoming 
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paper. 



6 Conclusions 

Solutions to the steady incompressible resistive magnet ohydro dynamics equa- 
tions in a curvilinear geometry are derived and discussed. These solutions 
describe a process where a sheared flow crosses a separatrix of the magnetic 
field and a curved current layer is formed in correspondence to the other sepa- 
ratrix. These solutions are the analogous in polar coordinates of the solutions 
in Cartesian coordinates discussed in [7] but with respect to the latter they 
present some new feature. In particular the configurations of the magnetic 
and velocity field in the curved geometry is much less symmetric than the 
one described by the solutions in Cartesian coordinates. Finally the possible 
applications of our solutions to the modeling of solar flares are mentioned. 
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